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Abstract 

We present a theoretical review of the possibilities offered by heavy flavour tagging in 
testing QCD features in 4-jet events at LEP1. Attention is focused mainly on the properties 
of the angular variables which are used in the experimental analyses to measure the colour 
factors of the strong interactions. We also comment on a possible settlement of the ongoing 
controversy about the existence in LEP1 data of SUSY events involving light gluinos. 
Integrated and differential rates of interest to phenomenological analyses are given and 
various tagging procedures are discussed. 
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1. Introduction 



The LEP1 era has now ended. Regarding QCD tests, the legacy passed on by the experiments 
ALEPH, DELPHI, L3 and OPAL is indeed a convincing indication that the SU(3) gauge group 
of the Standard Model (SAi) is the theory that describes the strong interactions By study- 
ing the process e + e~ — > Z — > hadrons decisive results have been achieved: the QCD coupling 
'constant' a s has been determined from jet rates and shape variables; both its flavour inde- 
pendence and running have been verified; multi-jet distributions reproduce QCD predictions to 
second order in a s ; many models proposed as alternatives to ordinary QCD have been ruled 
out; the triple gluon vertex coupling has been verified to be in agreement with QCD; finally, 
the colour factors that determine the gauge group responsible for strong interactions have been 
measured |2j. However, as far it concerns the latter measurements, somewhat less conclusive 
results have been obtained. 

In the 'colour factor analyses' in 4-jet samples, the basic idea is to measure the fundamental 
colour factors of QCD, that is, Ca and Cp and T F , which are determined from the SU(N C ) 
generators {T a )ij and from the structure constants f abc . Explicitly, one obtains Ca = Nc, 
C F = (N G - 1)/2N C , and T F = 1/2. In ordinary QCD (i.e., N c = 3): C A = 3 and C F = 4/3. 
In order to determine Ca, Cp and TpQ, one compares the theoretical predictions to the data by 
leaving the colour factors as free parameters to be determined by a fit. The comparison is made 
by resorting to differential spectra of variables based on angular correlations between jets, which 
are usually ordered in energy and where the two most energetic jets are 'identified' with the 
primary quarks (i.e., from the Z decay). Some of the most widely used are the Bengtsson-Zerwas 
angle Xbz, the Korner-Schierholz-Willrodt angle §ksw, the 'modified' Nachtmann-Reiter angle 
6* NR |§ and the angle between the two least energetic jets 6*34 (i.e., from the gluon splitting). 
The distributions in these quantities are different for QQgg and QQqq events. This is decisive 
as, in general, the main difference between QCD and alternative models is the predicted relative 
contributions of the two subprocessesQ This leads to different shapes of the angular distributions 
of the full 4-jet sample and thus to different predictions for the values of Ca, Cp and/or Tp 
(see, e.g., Ref. M). Ultimately, in fact, Cp, Ca and T F represent nothing else than the relative 
strengths of the interactions q — > qg, g — > gg and g — > qq, respectively. 

The present status of the colour factor analyses f|, || is that, although the experimental 
measurements agree well with ordinary QCD, it is not possible to rule out its Supersymmetric 
version. In fact, such a theory includes in the spectrum of the fundamental particles the spin 
1/2 partner of the gluon, the gluino g, which could be produced at LEP1 in QQgg events (via 
g* — > gg), provided it is not too heavy. Since gluinos are coloured fermions, their contribution 
would enhance the part of the 4-jet cross section with angular dependence similar to that of 
2Q2q events. Naively, one could alternatively say that N F is (apparently) increased, such that 
an experimental measurement could reveal SUSY signals in the form of an enhancement of Tr. 
In detail, gluinos with a mass ~ 2 GeV yield an expectation value for TpjCp within one standard 
deviation of the measured one ||. 

Like the QCD gauge boson, the gluino is neutral, it is its own anti-particle (i.e., it is a 
Majorana fermion) and its coupling to ordinary matter is precisely determined in terms of the 
QCD colour matrices and a s ||. However, whereas m g = 0, the gluino mass m g is an arbitrary 
parameter, and so is the lifetime Tg. In particular, in the theory it is natural for gluinos to be 
much lighter (a few GeV) than squarks if their mass is induced radiatively and if dimension-3 

1 Better, ratios among them: e.g., Ca/Cf and Tr/Cf, where Tr = NfTf, with Nf number of active flavours. 
2 E.g., the 2Q2q contribution to the total 4-jet rate at the Z peak is only about 5% in QCD but up to 30% in 
certain Abelian models 



SUSY breaking operators are absent from the low energy dynamics M. Thus, if gluinos are 
so light, they should be produced at LEP1 ||. Recently, there has been a renewed interest 
in this possibility, especially motivated by the (small) discrepancy between the value of a s 
determined by low energy deep-inelastic lepton-nucleon scattering and the one measured by 
the e + e~ CERN experiments 0. In this respect, it has been speculated that the evolution of 
a s is being slowed down by a contribution to the (3 function due to light gluinos. The latest 
experimental constraints still allow for the existence of relatively long lived and light gluinos, in 
the parameter regions pi: (i) m» ~ 1.5 GeV and rs ~ 10~ 8 s; (ii) rrig ~ 4 GeV and r„ ~ 10~ 10 s. 



In our opinion, if one adopts energy ordering of the jets (as it has been done so far), the 
effectiveness of the experimental analyses in putting stringent bounds on the measured values 
of Ca, Cf and Tr is largely reduced. In fact, the angular variables introduced above are most 
useful for emphasizing QCD features if one can distinguish between quark and gluon jets and 
assign the momenta of the final states to the corresponding particles^. It was this consideration 
that motivated two of us to advocate in Ref. |TT] the use of 4-jet samples with two of the jets 
clearly identified as originating from heavy flavour quarks. In this way, it is at least possible 
to distinguish some quark jets, namely those which originate from 6-quarks, from gluon jets, 
thus identifying the particles in the final states and defining the angular variables in the proper 
manner. A direct consequence of this approach is that one gains a greater discriminating power 
between events with different angular behaviour. Furthermore, an intrinsic advantage of the 
heavy flavour selection is that the QQgg event rates are reduced more (by a factor of 5) than 
the QQqq ones (which are divided by 3), such that their relative differences can be more easily 
studied. The final aim is to assess whether with 6-tagging techniques one is able to reduce the 
confidence level regions around the experimental measurements, in such a way to eventually 
shed new light on the possibility that SUSY signals are present in LEP1 data. 

Such techniques have been rapidly developed by the LEP1 collaborations in the past years 



12], and have become an excellent instrument to study heavy flavour physics. Their main 



features are well summarized, e.g., in Refs. [ 131 : in particular, we remind the reader that the 



method offering the best performances is the lifetime tagging by detecting a secondary decay 
/i-vertex. Yet, the crucial point in 6-tagging-based analyses is that the total statistics in 262jet 
events is greatly reduced compared to that available when no quark identification is exploited, 
not merely because only flavour combinations involving 6's are retained (and these are even mass 
suppressed |Hfl ) but also because the (single) tagging efficiency at LEP1 is rather low, e ~ 40% 
(for purities close to 1). Typically, one has to consider that the initial rate of 262jet events from 
pure QCD is of the order 0(8000) events per experiment: for pa 4 x 10 6 hadronic Z decays, 
T&b/rhad pa 22%, Ri pa 10% and e pa 30%, where is the 4-jet fraction, and considering two 
tags. Note that the statistics could well be increased by considering the data of all experiments 
altogether and/or by relaxing the requirements of high 6-purity (see discussion in Ref. ||1 1|| ). 
Indeed, experimental studies of 6-tagged 4-jet events have been already reported (showing 
that a considerable reduction of the systematics can be achieved) and others are in progress [JIB] . 

Independently of the results of the new analyses, there are however other possibilities offered 
by the /i- vertex devices to settle down the debate about colour factor measurements at LEP1. 
This is clear if one considers that long-lived gluinos might yield 4-jet events with detectable 
displaced vertices, as already recognised in Ref. JT7| and exploited in Ref. []TS[] . In the rest of 
this contribution, we briefly summarise the computational techniques, the tagging procedures 
and the main results already reported in Refs. jTT[ |TU| . 

3 When ordering in energy, one has to bear in mind that in Z — > QQgg decays for only half of the events the 
two lowest energy partons are both gluons |Q ! 



2. Calculation 



In carrying out the analyses described here we made use of the FORTRAN matrix elements already 
discussed in Ref. ]14| and presently used for experimental simulations ||, upgraded with the 



inclusion of the gluino production and decay mechanisms (see also Ref. [20]). The programs do 
not contain any approximations, the intermediate states 7* and Z being both inserted, and the 
masses and polarisations of all particles in the final states (of the two-to-four body processes) 
retained. The availability of the last two options is especially important if one considers, on the 
one hand, that in 6-tagged samples all final states are massive, and, on the other hand, that in 
proceeding to experimental fits one could well select restricted regions of the differential spectra 
of the angular variables, where the rates are likely to strongly depend on the spin state of the 
partons^]. 

In the analyses we are reporting about we have adopted four different jet-finding procedures: 
the JADE scheme (J) and its "E" variation (E), and the most recent Durham (D) and Geneva 
(G) algorithms ||21|| . We will present rates only for some of these, as none of the results drastically 
depends on the choice of the jet recombination scheme and/or the value of the jet resolution 
parameter, y cnt . As it concerns the exact definitions of the angular variables used here, these 
can be found in Ref. ||1 l||f|. Anyhow, to make clear the rest of the paper, we point out the 
following (with reference to the formulae of Ref. |jll|). When heavy flavour identification is 
implied, labels 1 & 2 refer to the two tagged jets whereas 3 & 4 to the two remaining ones. If 
no vertex tagging is assumed, jets are labelled according to their energy, E\ ^ E2 ^ E3 ^ E4.. 
The numerical values adopted for quark masses and SM. parameters can be found in Ref. [IPI 



3. Results 

3.1 Angular variables 

The first two figures illustrate the improvement that can be achieved by using heavy flavour 
tagging (two lower plots), compared to the results obtained when energy ordering is used (two 
upper plots). Distributions are normalised to one. The label 'Abelian' for the 2Q2g contribution 
refers to rates obtained in a particular Abelian model of QCD0, in which the self-interaction of 



gluons is absent [|23[ (see Ref. [|TTJ for details). Note that the QQqq subprocess is identical in 
both models. The following aspects must be emphasised in Figs. 1-2. 

(i) Concerning the angles xbz and ^*k$Wi the discrimination power between the behaviour of 
the QQgg and QQqq components slightly improves (more for the second variable). In the case 
of the first one, the original symmetric behaviour implied by the definition is restored. 

(ii) Another remarkable feature of Fig. 1 is that when 6-tagging is exploited, the differences 
between the Abelian and the non-Abelian behaviour of the 2Q2g component are dramatically 
emphasised, for both angles. In the case of the xbz variable these differences disappear when 
energy ordering is used. It should also be stressed that this is not the case for the $* KS w an gle 
only because we are using here an 'improved' definition of this variable (which shows sensitivity 

4 For this reason we have not used the results published in literature for the gluino decay rates, as these are 
averaged over the helicities of the unstable particle. Instead we have recomputed the relevant Feynman decay 
amplitudes by preserving the gluino polarisation and by matching the latter with the corresponding one in the 
production process Q. 

5 We stress that the &* KSW angle is actually a modification of the original one $ksw , proposed by two of us 
(see Ref. jll| for further details). 

6 We are aware that this model has been already ruled out in other contexts ]23] , but we regard it as a useful 
tool to demonstrate the sensitivity of the introduced angular variables to the features of QCD. 



to the gauge structure even without any flavour tagging), the original one being less sensitive 
to such differences. 

(hi) Regarding the 9* NR angle, once again the original symmetry is restored, such that this clearly 
helps to differentiate the 2Q2g and 2Q2q subprocesses. In contrast, the behaviour of the two 
different gauge interactions in the 2Q2g subprocess remains identical. 

(iv) Finally, for the angle #34, the improvement concerns the enhanced separation between the 
two different gluon behaviours, more than a discrimination between gluon and quark properties. 



3.2 Mass effects 

Fig. 3 shows the effect on the relevant differential spectra of the inclusion of the 6-quark mass 
in the computations. The angular distributions corresponding to 6-contributions are compared 
to those involving massless <i-quarks (q represents here a massless flavour). (Differences are 
even more dramatic when comparing bbbb and dddd distributions, though the former flavour 
combination corresponds to a suppressed contribution.) In particular, it is worth noticing that 
the (QCD) QQgg spectra is less affected by mass effects than the QQqq ones and that the 
Bengtsson-Zerwas and the (modified) Korner-Schierholz-Willrodt angles are more mass sensi- 
tive than the (modified) Nachtmann-Reiter and the angle between the two untagged jets. In 
general, it must be stressed that mass effects are not negligible and should be included in phe- 
nomeno logical analyses exploiting heavy flavour tagging. In contrast, when energy ordering is 
adopted (that is, when all quark flavours are retained in the four-jet sample), the above effects 
are largely reduced (see Ref. |14|)- 111 this case, in fact, the most part of quark combinations in 



the final state involves massless particles. Results have been shown for the case of the E scheme, 
as this is the algorithm for which mass suppression is somehow larger. However, similar effects 
are visible also in the other three cases: J, D and G. 



3.3 Gluino tagging 

When dealing with tagging a secondary vertex possibly due to gluino decays, several points 
must be addressed |19| . First of all, one has to confine oneself to secondary vertex analyses 



only[| however this technique has a larger efficiency than any other method [fl3| | . Second, the 
vertex has to be inside the detectors, so that only gluinos with t„ ~ 10~ 9 s can be searched 



for [0. Nonetheless, this represents an appealing opportunity, as a substantial part of the 
mentioned (m g , Tg) window could be covered by the experiment. In this respect, we exploit 
a sort of 'degeneracy' in lifetime between 6-quarks and gluinos, assuming that when making 
secondary vertex tagging one naturally includes in the 262jet sample also SUSY events, in which 
a g behaves as a b. We call such approach 'minimal trigger' procedure, as we propose a tagging 
strategy that does not take into account any of the possible differences between gluinos and 
6-quarks in 4-jet events with two secondary vertices (thus, in the following, we will generally 
speak of 'vertex tagging'). There are in fact at least three obvious dissimilarities, 
(i) Their charge is different, such that one could ask that the jet showing a displaced vertex 
has a null charge. This would allow one to isolate a sample of pure SUSY events. However, we 
remind the reader that to measure the charge of a low energy jet in 4-jet events has not been 
attempted before and that the procedure would certainly have very low efficiency (in isolating 
a very broad hadronic system in an environment with high hadronic multiplicity). 

7 Thus neglecting other forms of heavy flavour tagging: such as the high pt lepton method, as gluinos do not 
decay semileptonically. 



(ii) Gluino lifetimes much longer than that of the b are still consistent with experiment (note 
that r b » 10 -12 s), such that recognising a displaced vertex with decay length d 3> 3 mm (that 
of the b) would allow one to immediately identify gluinos. Unfortunately, most of the 6-tagged 
hadronic sample at LEP1 has been collected via a bi-dimensional tagging (see, e.g., Ref. ||24j| ). 
Thus, different d's could well appear the same on the event plane. Furthermore, tagging a d > 3 
mm vertex would allow one to separate gluinos with Tg > r b , but this would not be helpful if 

Tg < T b . 

(hi) Other than in lifetime, 6-quarks and gluinos can differ in mass as well, such that one might 
attempt to exploit mass constraints to separate SUSY and pure QCD events. However, on the 
one hand, one could face a region of m b -m g degeneracy and, on the other hand, one should cope 
with the ambiguities related to the concept of mass as defined at parton level and as measured 
at hadron level. 

We emphasise that measuring the charge of the vertex tagged jet, attempting to disentangle 
different decay lengths or measuring partonic masses could well be further refinements of the 
procedure we are proposing, and that their exploitation certainly does not spoil the validity of 
the latter. In addition, all these aspects necessarily need a proper experimental analysis, which 
is beyond the scope of a theoretical study. 

The steps of our analysis are very simple. Under the assumption that fe's and <fs are not 
distinguishable by vertex tagging, one naturally retains in the sample all SUSY events (whereas 
the ordinary QCD components are reduced as mentioned earlier). Then, it is trivial to notice 
that there exist clear kinematic differences between the QQgg, QQqq and QQgg components. 
This is shown in Fig. 4, where we plot the quantities = Mfj/s = (pi + Pj) 2 /s, where ij = 12 
or 34 and s = E^ m = {M z ) 2 ■ The behaviour of the curves is dictated by the fact that gluinos 
are always secondary products in 4-jet events, whereas quarks and gluons are not (lower plots). 
When no vertex tagging is exploited, such differences are washed away (upper plots). The value 
irtg = 5 GeV is assumed for reference, the shape of the distributions being qualitatively the 
same regardless of it. Therefore, if one simply asks to reject events for which, e.g., Y i2 > 0.2 
and/or Y" 34 < 0.1, one gets for the total rates of the three components the pattern exemplified 
in Fig. 5. From the drastic predominance of 2Q2g events in the complete 'unflavoured' sample 
(no vertex tagging, top left), one first obtains that the total rates of ordinary QCD events are 
significantly reduced compared to those of SUSY events (after vertex tagging, top right), and 
eventually that the QQgg fraction is always comparable to that of QQgg + QQqq events (when 
also the kinematic cuts are implemented, bottom left): and this is true independently of the 
gluino mass, of the jet algorithm and of the y cut value used in the analysis (see Ref. jHfl). 

Therefore, after our event selection, SUSY signals would certainly be identifiable, as a clear 
excess in the total number of 4-jet events with two displaced vertices. Thus, the final goal will 
be to definitely assess the presence of gluinos at LEP1 or contradict the latter, at least over 
appropriate regions of masses and lifetimes^. We finally stress that, as we are concerned here 
with total rates and not with differential distributions, the event number should be sufficient 
to render the analysis statistically significant^ and that the usual ambiguities related to the 
fact that gluino effects on the total number of 4-jet events are comparable in percentage to the 
systematic uncertainties due to the jet hadronisation process and/or the y CVLt selection procedure 

8 And this should certainly be done after the appropriate MC simulations, including the details of the detectors 
and of the tagging procedure as well as a generator where ra g and Tg enter as free parameters to be determined 
by a fit. 

9 In this respect we acknowledge that many of the aspects of our approach were already employed in 
Ref. Q, the tagging procedure sketched there being however well beyond the statistical possibilities of the 
LEP1 experiments. 



are much less severe in our approach^ However, since the key point of the present approach is 
to exploit the b/g vertex degeneracy, a highly enriched heavy flavour sample should be selected 
in this case. 

3.4 Gluino decays 

Before closing, we should mention that a further aspect must be kept into account when at- 
tempting our analysis. It concerns the kinematics of the gluino decays. In the most sponsored 
SUSY framework ||, the dominant gluino decay modes are g — > qq^ and g — > (77, where 7 
represents a 'photino' (better, the Lightest Supersymmetric Particle, which is a superposition 
of the SUSY partners of the neutral gauge bosons of the theory). Furthermore, the qcpi channel 
is, in general, largely dominant over the #7 mode 0. 

The crucial point is that in both cases the gluino decays into a jet with missing energy. It 
is not our intention to discuss the possibility of selecting such a signature, as we are mainly 
concerned here with the fact that the energy left to the hadronic system Eh is above the ex- 
perimental cuts in minimal hadronic energy, which are used to reduce the backgrounds (e.g., 
in Ref. || the threshold was set equal to 3 GeV). In Fig. 6 (first three plots) we show the E h 
spectra after the gluino decay, in both the channels. Two kinematic decay configurations are 
considered: a massless photino, and a massive one (i.e., = \j1m-g). The message is that in 
the most likely SUSY scenario (i.e., three-body decay dominant and massless photino) all gluino 
events should be retained in the event selection. Conversely, Fig. 6 illustrates the percentage 
of these which will pass the adopted trigger requirements. Finally, in the bottom right plot of 
Fig. 6 we show the dependence of the SUSY rates on the value of trig. Below m g « 5 GeV, the 
mass suppression is always less than a factor of 2. 

4. Summary and conclusions 

In this paper we have stressed the importance of using at LEP1 samples of 4-jet events, in which 
two of the jets show a displaced vertex. This could help to improve some of the experimental 
tests of QCD and to possibly settle down the ongoing dispute about the existence of SUSY 
events in the data. Those presented are theoretical results, which should be in the end verified 
by detailed MC simulations. However, it is our opinion that they indicate that the matter raised 
and procedures similar to the ones outlined here would deserve experimental attention. 

We are grateful to Ben Bullock for reading the manuscript. We also acknowledge the use of 
a FORTRAN code produced by Kosuke Odagiri to evaluate the gluino decay rates. This work is 
supported by the MURST, the UK PPARC, the Spanish CICYT project AEN 94-0936, the EC 
Programme HCM, Network "Physics at High Energy Colliders", contracts CHRX-CT93-0357 
DG 12 COMA (SM) and ERBCHBICHT (RMT), and by the Graduiertenkolleg "Teilchen- 
physik" (JBT). 



References 

[1] See, e.g.: T. Hebbeker, Phys. Rep. 217 (1992) 69, and references therein. For later refer- 
ences, see citations made in Ref. ||TT] . 



10 These are in fact the underlying difficulties of any analysis based on the 'unflavoured' hadronic sample and/or 
the jet energy ordering, which have not been overcome even in recent improved studies p5[. 



[2] ALEPH Coll., Phys. Lett. B284 (1992) 151; DELPHI Coll., Phys. Lett. B255 (1991) 466; 
DELPHI Coll., Z. Phys. C59 (1993) 357. 

[3] M. Bengtsson and P.M. Zerwas, Phys. Lett. B208 (1988) 306; J.G. Korner, G. Schierholz 
and J. Willrodt, Nucl. Phys. B185 (1981) 365; M. Bengtsson, Z. Phys. C42 (1989) 75; O. 
Nachtmann and A. Reiter, Z. Phys. C16 (1982) 45. 

[4] S. Bethke, A. Ricker and P.M. Zerwas, Z. Phys. C49 (1991) 59. 

[5] OPAL Coll., Z. Phys. C65 (1995) 367. 

[6] For review papers, see e.g.: H.P. Nilles, Phys. Rep. 110 (1984) 1; H.E. Haber and G.L. 
Kane, Phys. Rep. 117 (1985) 75; R. Barbieri, Riv. Nuo. Cim. Vol. 11 No. 4 (1988) 1. 

[7] G.R. Farrar and A. Masiero, preprint RU-94-38, October 1994; G.R. Farrar, preprint RU- 
94-70, August 1994. 

[8] P. Nelson and P. Osland, Phys. Lett. 115B (1982) 407; B.A. Campbell, J. Ellis and S. 
Rudaz, Nucl. Phys. B198 (1982) 1; G.L. Kane and W.B. Rolnick, Nucl. Phys. B217 (1983) 
117; B.A. Campbell, J.A. Scott and M.K. Sundaresan, Phys. Lett. 126B (1983) 376. 

[9] L. Clavelli, Phys. Rev. D46 (1992) 2112; L. Clavelli, P. Coulter and K. Yuan, Phys. Rev. 
D47 (1993) 1973; M. Jezabek and J.H. Kiihn, Phys. Lett. B301 (1993) 121; R.G. Roberts 
and W.J. Stirling, Phys. Lett. B313 (1993) 453; M. Carena, L. Clavelli, D. Matalliotakis, 
H.P. Nilles and C.E.M. Wagner, Phys. Lett. B317 (1993) 346; T. Hebbeker, Z. Phys. C60 

(1993) 63; D.V. Shirkov and S.V. Mikhailov, Z. Phys. C63 (1994) 463; J.L. Lopez, D.V. 
Nanopoulos and X. Wang, Phys. Lett. B313 (1993) 241. 

[10] G.R. Farrar, Phys. Rev. D51 (1995) 3904; B. Kileng and P. Osland, Z. Phys. C66 (1995) 
503. 

[11] S. Moretti and J.B. Tausk, Z. Phys. C69 (1996) 635. 
[12] For a list of relevant references, see Ref. [|Il[ . 

[13] S. Squarcia, preprint CERN-PPE-94-69, April 1994; A. De Angelis, preprint UDPHIR 
95/01/AA, January 1995; Z. Fodor, Phys. Rev. D41 (1990) 1726, and references therein. 

[14] A. Ballestrero, E. Maina and S. Moretti, Phys. Lett. B294 (1992) 425; Nucl. Phys. B415 

(1994) 265; Proc. of the XXIXth Rencontres de Moriond, Meribel, France, March 1994. 

[15] K. Hamacher, talk given at the XXVIIIth International Conference oh High Energy Physics, 
25-31 July 1996, Warsaw, Poland, to appear in the Proceedings. 

[16] S. Bentvelsen, private communication. 

[17] J. Ellis, D.V. Nanopoulos, and D.A. Ross, Phys. Lett. B305 (1993) 375. 
[18] F. Cuypers, Phys. Rev. D49 (1994) 3075. 

[19] S. Moretti, R. Munoz-Tapia and K. Odagiri, preprint UGR-FT-67/96, Cavendish-HEP- 
96/6, DFTT 25/96, August 1996 (in preparation). 



[20] R. Munoz-Tapia and W.J. Stirling, Phys. Rev. D49 (1994) 3736. 



[21] JADE Coll., Z. Phys. C33 (1986) 23; Phys. Lett. B213 (1988) 235; N. Brown and W.J. 
Stirling, Phys. Lett. B252 (1990) 657; Z. Phys. C53 (1992) 629; S. Bethke, Z. Kunszt, D.E. 
Soper and W.J. Stirling, Nucl. Phys. B370 (1992) 310. 

[22] OPAL Coll., Phys. Lett. B235 (1990) 389; S. Bethke, in Proc. of the XXVth Rencontres 
de Moriond, Les Arcs, Savoie, France, March 1990. 

[23] K.J.F. Gaemers and J. A.M. Vermaseren, Z. Phys. C7 (1980) 81; O. Nachtmann and A. Re- 
iter, Z. Phys. C16 (1982) 45. 

[24] See, e.g.: OPAL Coll., Nucl. Instr. Meth. A348 (1994) 409. 



[25] A. de Gouvea and H. Murayama, preprint LBL-39030, June 1996. 



e + e -> j j j j 



0.011 



y V = 0.0015 

q C B i <QCD) cut 
q q g g {Abelian) 



0.009 - 



0,007 - 



0.005 - 



m 
4) 

Id 0.003 
J3 




cs 
m 

X 

^ 0.009 

b 

T3 



0.007 - 



0.005 



0.003 p 



0.001 



I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 



Q = q 



i I ■ I i I ■ I ■ I i I ■ I i I i 



Q = b 




i ■ i ■ i ■ i ■ i i i ■ i i i 



SQQ qq 

I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 



Q = q 



■ I ■ I ■ I ■ I ■ I ■ I i I ■ I i 



Q = b 




Q 

■e- 

MS „ 
* * 



ft 

CD 
ft 



20 40 GO 60 100120140160 20 40 60 80 100 120 140 160 180 



Xbi (degrees) 



KSW 



(degrees) 



Fig. 1 



Distributions in the angles \BZ and &*ksw f° r the QQgg (QCD and Abelian) and QQqq components 
in 4-jet events, in the D scheme with y^ ut = 0.0015, without (Q = q) and with (Q = b) 6-tagging. 
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Fig. 2 



Distributions in the cosine of the angles 9* NR and #34 for the QQgg (QCD and Abelian) and QQqq 
components in 4-jet events, in the D scheme with y® ut = 0.0015, without (Q = q) and with (Q = b) 
6-tagging. 
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Fig. 3 



Distributions in the angles xbz and &ksw anc ^ m ^ ne cos i ne °f the angles 9* NR and #34, for the QQgg 
and QQqq components in 4-jet events, in the E scheme with yf ut = 0.01, where Q = b,d and q is a 
massless quark. 
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Distributions in the rescaled invariant masses Yij = M^/s, where ij = 12,34, for ordinary QCD and 
for SUSY 4-jet events, in the D scheme with y^ t = 0.002, without (Q = q) and with (Q = b) vertex 
tagging. Here, nig = 5 GeV. 
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Fig. 5 



Total cross sections of ordinary QCD and of SUSY 4-jet events, in the D scheme, without (Q = q) and 
with (Q = b) vertex tagging, and after the kinematic cuts, for three different values of m~ g . 




Fig. 6 



Differential distributions in the hadronic energy of the 'gluino jet' after the two possible SUSY decays, 
in the D scheme, for various combinations of m g and m^; and total cross section of gluino events in 
4-jets, as a function of m g and for three different values of y^ut- 



